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ABSTRACT 

Recent progress in observational studies of magnetic activity in M dwarfs urgently requires 
support from ideas of stellar dynamo theory. We propose a strategy to connect observational 
and theoretical studies. In particular, we suggest four magnetic configurations that appear rel¬ 
evant to dwarfs from the viewpoint of the most conservative version of dynamo theory, and 
discuss observational tests to identify the configurations observationally. As expected, any 
such identification contains substantial uncertainties. However the situation in general looks 
less pessimistic than might be expected. Several identifications between the phenomenology 
of individual stars and dynamo models are suggested. Remarkably, all models discussed pre¬ 
dict substantial surface magnetic activity at rather high stellar latitudes. This prediction looks 
unexpected from the viewpoint of our experience observing the Sun (which of course differs in 
some fundamental ways from these late-type dwarfs). We stress that a fuller understanding of 
the topic requires a long-term (at least 15 years) monitoring of M dwarfs by Zeeman-Doppler 
imaging. 
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1 INTRODUCTION 


The intensive investigation of stellar magnetic activity has already a 
history of about 40 years, and many important results have been ob- 


tained 

Hartmann & Noyes 1987 

Berdvugina 2009; Saar 20Ut 

Donati 

2011: Vidotto etaU2014 

; Kovari & Olah 2013). On one 


hand, this activity is believed to be similar to some extent to the so¬ 
lar cyclic activity. In particular, the magnetic fields of many types 
of stars are believed to be excited by stellar dynamo action, driven 
by the joint action of differential rotation and an additional mirror- 
asymmetric driver which produces poloidal magnetic field from 
toroidal. The physical nature of this mechanism is still debated for 
the Sun, and there are even more uncertainties for other stars. How¬ 
ever we have no intention to enter this discussion now and for the 
sake of simplicity refer to this factor as the ff-effect, whatever its 
physical nature. 

On the other hand, magnetic activity in any particular star can 
look quite distinctively different from the solar case: starspots on 
the most investigated magnetically active stars are much larger than 
on the Sun, covering a larger part of the surface than sunspots, etc. 
These specific features should somehow be related to the properties 
of stellar dynamos. However in practice any such relation is ill- 
understood. 

In particular, the relation between the magnetic phenomenol¬ 
ogy of M dwarfs and the underlying stellar dynamo action is a topic 
of continuing discussions at scientific meetings. Experts in obser¬ 


vational investigations of starspots on M dwarfs have stressed many 
times that some input from dynamo theory to the problem is very 
desirable. There is a general feeling that there is a basic difference 
between dynamo action in M dwarfs which are fully convective 
stars, and in the Sun which has a relatively thin convective shell. 
As far as we can recall, it was E. Ergma who attracted the attention 
of one of us (Dmitry Sokoloff) to the problem as early as 1993. 


There are many theoretical models that explore various de¬ 
tails of stellar dyn amo action for a range of stellar types, in clud- 
ing M dwarfs (e.g. Jchabrier & K ii ke r 1 200fil : 1 B r() w n i rig 1 2008h . 3D 
dynamo models give very important information concerning mag¬ 
netic activity of fully convective stars. In particular, the models 
clarify spot formation mechanism jYadav etaL 2015^ a nd pr edict 
the location of dipole-dominated M-dwarfs (Schrinne r et al. .120121: 


iGastine et al. 11200 , 120131 : ISchrinner et al. Il2014ll in terms of their 

rotation velocities and convection efficiency, giving in this way an 
observational test for verification of the model. However, some¬ 
times the predicted magnetic configurations depend significantly on 
various small and quite uncertain factors, such as the boundary con¬ 
ditions at the stellar surface; it is thus very difficult to adequately 
reproduce all these details in a particular dynamo model. 


Therefore, although considerable progress has been made in 
developing comprehensive 3D MHD simulation s of fully or nearly 


fully convective stars (e.g. Dobler et al. 2006j; Browning 

2008; 

Schrinner et al. 2012; Gastine et al. 20121; Schrinner et al. 

20ll; 

Yadav et al. 2015!), these models are not yet fully satisfactory. For 
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instance, applied to the Sun the 3D models still have difficulties in 
tuning parameters to reproduce the observed butterfly diagram. 

Depending on assumptions and approximations used in nu¬ 
merical simulations, modern dynamo models may predict dif¬ 
ferent dymano states. For instance, for fully convective stars it 
is often suggested that dynamos will be of a 2 type, and of¬ 
ten non-axisymmetry is predicted. However, the steady a 2 solu¬ 
tions seem to be not fully in accordance with observations, where 
axisymmetric and non-axisymmetric configurations are found in 
very similar objects (see Sect. [4] also lGastine et al. 1 120141) . Note 
that relatively weak toroidal field is often taken as an indicator 
of non- axisy mme try. O n the other hand, as was shown in e.g., 
[Schrinner et al. 1 I2012} ). oscillating or-dynamos are also possible 
under certain conditions. 

Furthermore, very recent dynamo models predict that dynamo 
can be bistable, the eventual stable state depending on initial condi¬ 
tions. This could explain the observed distributions of the magnetic 
field geome tries among objects with similar parameters (see e.g., 
ISchrinner et al. 120141) . 

One of the difficulty in using 3D models for the extensive ex¬ 
ploration of stellar dynamos is the requirenment of large computa¬ 
tional resources, making it difficult to use outside of a substantial 
programme, and limiting possible exploration of parameter space. 
As alternative, the “hybrid” models that use transport coefficients 
deduced from MHD simulations in m ean field models may provide 
a more satisfactory way for ward fe.g jDube & Charbonneaull 20131 : 
IPipin & Kosovichev 120131 1. Nevertheless, in spite of a wide variety 
of available models, none of them provides a complete understand¬ 
ing of dynamo action in low mass stars. 

The above considerations suggest that progress in theoretical 
modelling and observations has led to a rather strange situation 
where quite amazing advances in the last 40 years have brought 
rather less physical understanding of the problem than might have 
been anticipated. Our feeling is that coordination of efforts be¬ 
tween observers and dynamo modellers is urgently needed to fur¬ 
ther progress in the topic. 

We envisage the format of the coordination to be as follows. A 
risk has to be taken and a particular approach to dynamo modelling 
has to be accepted as a starting point. A list of magnetic configu¬ 
rations that can be expected for M dwarfs in the framework of this 
basic model has to be constructed, together with observational tests 
capable of distinguishing the configurations observationally. If pos¬ 
sible, the existing observational data have to be used to identify the 
configurations, and the perspective of future observations should 
be discussed. The point is that the full inspection of the list may 
(and will) be quite demanding on telescope time and other mate¬ 
rial resources, and the community will have to make a reasonable 
decision on how to manage resources in order to achieve this goal. 
Further steps in development of the model would have to be moti¬ 
vated by particular difficulties in interpretation of observations. 

The aim of this work is to explore the magnetic configura¬ 
tions that can be excited by a more or less standard mean-field dy¬ 
namo model based on differential rotation (D-effect) and mirror- 
asymmetric convective motions (the a-effect), and then to follow 
the above plan as far as it is possible at the moment. We appre¬ 
ciate that several important points of our (as any other) dynamo 
model are debatable and can only, at best, partially represent the 
true physics. In particular, we recognize that in other models the 
differential rotation may be substantially reduced compared to the 
model we consider in this work. We stress however, that this pa¬ 
per is more a demonstration of principle than an attempt to of¬ 
fer a definitive resolution of any particular issue, and should be 


viewed as such. We do not attempt to give a complete review of 
the field, but rather just to indicate some key points and references. 
Of course, the parameters and details of the models are likely to be 
modified in course of progress in understanding of the topic. 

This paper is in some sense related to the recent paper by 
I Kitchatinov et al. I d2014l ) who, however, approached the topic from 
a quite different viewpoint. They suggested an explanation for the 
observed phenomenology of magnetic activity of M dwarfs in a 
way analogous to the solar case: that explanation does not require 
the idea of dynamo bistability. Here, however, we do not exclude 
the latter in principle and address the problem from an alterna¬ 
tive viewpoint, and discuss the richness of dynamo models for M 
dwarfs that have more or less comparable complexity. As we will 
show below, four rather than two configurations look possible and 
we suggest how to clarify which (if any) are really of potential in¬ 
terest, and whether something more complex is needed. 


2 THE DYNAMO MODEL 


Dynamo theory explains stellar magnetic activity by two basic ef¬ 
fects: i) poloidal field production from toroidal field by cyclonic 
motions (the a-effect) and ii) conversion of poloidal field back to 
toroidal by the action of differential rotation (the D-effect) and pos¬ 
sibly by the a-effect as well. The efficiency of the two effects in 
generating magnetic fields are measured by the dimensionless pa¬ 
rameters 


*7t 

for the D-effect (AD is a measure of the angular velocity variation 
within the convection zone and r/ T is the eddy magnetic diffusivity), 
and 


C a 


n T 


( 2 ) 


for the a-effect (Krause & Radler 1 980l) , where ao is a representa¬ 
tive value. 

For our illustrat i ve pur poses we use the same dynamo model 
as IlCitchatinov et al. 1 l2014h and therefore describe it only briefly. 
The model is of a 2 D-type, i.e. taking the toroidal field genera¬ 
tion by the a-effect into account. The differential rotation in M 
dwarfs is small feames et al. 112005 ) and it is often neglected in 
models of their dynamos, thus giving a 2 -dynamos. We note that 
it is possible to obtain similar magnetic field configurations with 
or without differential rotation. For instance, axisymmetric fields 
have been obtained in rece nt simulations with pure a 2 -dynamos 
(e.g.. ISchrinner et al, I 201_2|). On the othe r hand , mean-field mod¬ 
els of a 2 -dynamos by Chabrier & Kiikerl ({2006) (which allow for 
the rotationally induced anisotropy of the a-effect) predict non- 
axisymmetric global fields similar to an equatorial dipole. But if 
differential rotation is present, it may oppose the non-axisymmetry 
by converting azimuthal variations into small-scale radial or merid¬ 
ional structure, enhancing diffusive decay dMoss & Brandenburg I 
11995ll . Because the results of simulations are very sensitive to ap¬ 
proximations and assumptions made, in the mean-field models that 
we use in this paper we explicitly include the contribution from 
differential rotation. The point here is that although AD decreases 
with decreasing stellar mass, the eddy magnetic diffusion r] T de¬ 
creases even faster so that the parameter Cn (Eq.Q]l increases. This 
trend was_found in mean-field models of differential rotation by 
iKitchatinov & Olemskov IJ20I ill . Note again that this model does 
not prescribe the eddy transport coefficients but estimates them in 








































terms of the entropy gradient, the entropy being a dependent vari¬ 
able of the model. In particular, the value v T ^ 1.2 x 10 11 cm 2 s _1 
was estimated for the middle radius r = R star /2 of a 0.3M 0 star ro¬ 
tating with a period of 10 days (v T varies moderately with depth). 
We shall use the differential rotation computed for this fully con¬ 
vective star in our dynamo model. The eddy magnetic diffusion is 
of the same order of magnitude as the eddy viscosity l l Yousef et al. I 
120031 ), so that the magnetic Prandtl number. 


Pm = v T /ty T 


(3) 


does not differ much from unity. Taking P m = 1 gives Cn - 290 
which is not small in spi te of a rather small diffe rential rotation of 
about 1.5% (see Fig.l in lKitchatinov et al. 1120141) . 

The value of C a (Eqn.0 is more difficult to estimate. The ori¬ 
gin of the a-effect is not certain for the Sun iCharbonneau II 20 I 4 I) 
and even less so for M-dwarfs. Order of magnitude estimation is 
nevertheless possible. The Rossby number Ro = ProjT is small 
for the 0.3M 0 star considered with P rot = 10 days (r is the con¬ 
vective turnover time with a characteristic value of a few months 
(e.g jBarnes & KirrfeOlOl )). The a 0 parameter for this case of con¬ 
vection strongly affected by rotation is expected to be close to its 
maximum possible value of the convective velocity u c ~ 2r\t£\ £ 
being the mixing length. The estimate C a ~ 3R staI /£ follows from 
Eq.0. With the plausible value of R star /£ ~ 10, this gives C a ~ 30, 
which is not much above the critical C a values for onset of dynamo 
action in our model (a stellar structure model gives R stsI /£ = 4.3 for 
the middle radius r = R staI /2 leading to a somewhat smaller C a ). 

It may be expected, however, that back reaction of gen¬ 
erated magnetic fields oiy motion^ affects primarily the a-effect 
dBrandenburg & Subramanian 1120051) . The majority of mean-field 
dynamo models employ the a-quenching as the only nonlinearity. 
The quenching reduces C„ to close to the threshold value for onset 
of dynamo action. 

Note that estimates of C a and Cn suggest the importance of 
differential rotaion only in the models that we have chosen to use 
in this paper. This, of course, does not exclude the pure a 2 -dynamo 
action in M dwarfs, as was stressed above. However, the theoretical 
ability to predict dynamo governing parameters does not allow one 
to obtain firm predictions and we believe that observations may 
play a decisive role here. 

We apply our kinematic dynamo model to estimate the thresh¬ 
old values of C a for magnetic fields of different equatorial and axial 
symmetries. Our model assumes uniform diffusion and that the a- 
effect is uniform with radius and varies as cos 9 with co-latitude 
9, a = ao cos 9. The model solves the ei genvalue problem for the 
mean-field dynamo equations (cf., e.g., [Krause & Radler 1119801 ) 
numerically by applying a grid point method in radius and a lat¬ 
itudinal expansion in Legendre polynomials. Computations were 
performed on a uniform radial grid of 301 points, and with 40 
Legendre polynomials in the latitudinal expansion. Vacuum con¬ 
ditions were applied at the outer boundary. These conditions de¬ 
mand that the toroidal field vanishes at the surface. We cannot 
exclude however the possibility that the vacuum boundary condi¬ 
tion is not the whole story. It is possible to have a field above the 
surfac e whic h is force-free, with finite toroidal field, e.g. iMestel I 
1 1966|), later Milsom & Wright I (1 19761) , see also the discussion in 
iMoss & Sokoloff I d2009h. ~ 

Our model takes a to be independent of radius. Arguing from 
a possible correlation with the kinetic helicity (which is often used 
as a proxy for the a-effect), it might be appropriate to take a as an 
increasing function of radius. As stressed above, we only wanted 
to provided an illustrative model on which to hang our arguments. 
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rather than to provide an extensive exploration of parameters, and 
so did not pursue this. 

For numerical reasons the model requires an inner radial 
boundary. This artificial boundary was imposed at r ; = 0. lR star 
(Rstar = 212 Mm). Conditions for an interface with a perfect con¬ 
ductor were imposed at the inner boundary. 

Our toy model has a distinctl y non-cylindrical rotation profile 
(Fig. 1 of iKitchatinov et al. I2OI4I) . in contrast to predictions of 3D 
simulations of r apidly rotating models with strong dynamo action 
(e.g. lBrowning 1120080 . Feedback from Lorentz torques can reduce 
differential rotation. Weaker fields will merel y modify the underly¬ 
ing rotation field fe.g. lMoss & Brooke1l2000h . 

We appreciate that a number of potentially significant 
changes connected with the inclusion of non-linear effects (alpha¬ 
quenching, Lorentz force, etc.) could be made to the model and 
this could affect the saturated states (e.g.. IMoss & Brooke 200(1 : 
[Browning l2008i) . To avoid these complications, below we consider 
only linear models. As was stressed above, we prefer to start our 
exploratory work with simple models that will be advanced once 
required by the analysis of observed data. 


3 DYNAMO EXCITED MAGNETIC FIELD 
CONFIGURATIONS 

We use the stand ard notations S m and A m for global modes of mag¬ 
netic fields (cf. lKrause & Radler II 1980h , where the first letter S or 
A denotes field configurations that are symmetric or antisymmetric, 
respectively, with respect to the equator. In this notation, m is the 
azimuthal wave number. A0 and SO, therefore denote axisymmetric 
global fields, which can be steady or oscillatory. Al and S1 are non- 
axisymmetric modes with m = 1; each field component changes 
sign twice over an entire longitudinal circle. The modes with in ^ 2 
normally have dynamo excitation thresholds that are too high for us 
to consider them as a probable outcome of global dynamo action 
in M dwarfs. All the global modes combine poloidal and toroidal 
field components. The meaning of toroidal and poloidal fields for 
axisymmetric modes is well known. Lines of toroidal fields in non- 
axisymmetric modes lay on concentric spherical surfaces. Poloidal 
non-axisymmetric fields have toroidal vec tor pot entials and are 
supported by toroidal currents (see lChandrasekhar 11 96(1 for more 
details). 

Thinking straightforwardly, we have to conclude that the 
modes that can be excited for the lowest C„, i.e. the Al mode for 
P m < Ll and SO for P m > 1.1, are the modes likely to be ex¬ 
cited by dynamo action in M- dwarfs. The point however is that 
the difference in excitation condition between the modes SI and 
Al for lower P m and between the modes A0 and SO for larger P m 
respectively is quite small and it would be unrealistic to think that 
contemporary dynamo models are sufficiently near physical reality 
to support such distinctions. Thus we allow that any one of these 
modes might be excited in a given M dwarf. We assume that the 
probabilities of getting any of the configurations for a given P m are 
comparable for the members of the pair. 

To be definite, we ignore at least for the time being the non¬ 
oscillating modes shown by dots in Fig. Q] We stress that the os¬ 
cillating nature of the nonaxisymmetric modes SI and Al is non¬ 
physical: there is a rotating frame in which each of these magnetic 
configuration is non-oscillating. In other words, oscillation of the 
modes SI and Al is just a result of longitudinal drift of the field 
patterns, viewed from a given line of sight. 

The period of activity cycles for truly oscillating modes are 
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Figure 1. A dynamo model for a M dwarf: marginal values of C a (Eqn.0 
for axisymmetric modes (AO) and non-axisymmetric m = 1 modes (Al) 
of dipolar symmetry, and axisymmetric modes (SO) and non-axisymmetric 
modes (SI) of quadrupolar symmetry. Dots indicate steady axisymmetric 
modes and the solid line shows the oscillatory modes. 



Figure 2. The period of activity cycles for truly oscillating modes (to be 
compared with the 22-year solar activity cycle). 


shown in Fig. [2] Butterfly diagrams for these modes are shown in 
Fig. [3] Note that the activity waves propagate polewards. Q 

As butterfly diagrams are irrelevant to nonaxisymmetric 
modes, we represent their field configurations by a diagram in 
latitude-longitude coordinates (Fig. 0 

Of course, such figures have previously been presented in the 
literature, but perhaps with not so much detailed inter-comparison. 
Our aim is to discuss to what extent the four dynamo modes under 
discussion can be distinguished observationally at present time or 
in the visible future and how it could in principle be done. 


1 Application of the same model to a sun-li ke star gave equatorward mi¬ 
gration (see Fig. 5 in Kitchatinov et al. ] ioooh . 
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Figure 3. Butterfly diagrams for truly oscillating modes: upper pair - mode 
SO, lower pair - mode A0. Patterns of radial field at the surface and toroidal 
field at small depth are shown. Pm = 2.4. 


4 DISTINGUISHING DYNAMO MODES 
OBSERVATIONALLY 

First of all, no single mode from the above list has a configuration 
similar to the solar case, i.e. an oscillating axisymmetric mode of 
dipolar symmetry with equatorward propagation, and with toroidal 
magnetic field (i.e. stellar spots) located near to the stellar equator. 

The following point looks instructive. The toroidal magnetic 
field becomes strong near to the stellar equator only for the mode 
SO. In other words, confirmed observation of stellar spots near the 
stellar equator is a diagnostic feature for the mode SO. Note how¬ 
ever that, as can be seen from Fig. 3, more or less strong equatorial 
spots can also be observed for A0 configurations from some incli¬ 
nations. The exact latitudinal pattern of toroidal field is likely to 
vary with the angular form of a, e.g. a oc sin 2 6/cos 8 would push 
it towards the equator (simply because the maximum of one of the 
dynamo drivers becomes closer to the stellar equator). In density- 
stratified convection, it seems likely that helicity (which is often 
used as a proxy_of the a-effect) is indeed stronger near to the equa¬ 
tor I lGastine et al. 1 120121) . Excitation of mixed parity solutions is 





























































Dynamos in M dwarfs 5 


RADIAL FIELD 



TOROIDAL FIELD LINES 




0 100 200 300 

LONGITUDE 



0 100 200 300 

LONGITUDE 


Figure 4. Configuration of nonaxisymmetric modes in latitude-longitude 
coordinates: upper pair - mode Al, lower pair - mode SI. Pm = 1.3. 


known for nonlinear spherical dynamos (e.g. iBrandenbur g et al, I 
1 19891 : Ijennings & Weissl 1 199 lh and for solar activity at th e end 


of the Maunder minimum I SokolofF & Nesme-Ribes 11 1994h . and 
would also probably complicate interpretation. 

Modes SI and Al have a pronounced longitudinal modula¬ 
tion of their nonaxisymmetric toroidal fields. We note that photo¬ 
metric data for the Ca lines obtained in the framework of the H- 
K project allow recogn ition of rotatio n modulation for some stars 
dBaliunas et al. I Il995l) . Remarkably dKatsova et al. I l2010h . this 
modulation is well-pronounced not for the stars where the activity 
cycle is strong (“excellent” in the classification of iBaliunas et al. I 
l ll995l) ). but rather for the stars where cycles are moderate (“good” 
of Baliunas et al.’s classification). If such modulation were to be 


isolated for flare and flat stars in the classification of lBaliunas et al. I 
d 19951) . i.e. for stars without marked cycles, these stars could be 
considered to have magnetic fields of SI or Al configuration. 

In order to distinguish between S and A modes it is neces¬ 
sary to compare the magnetic field direction in “Northern” and 
“Southern” hemispheres. The problem is that to perform this test 
observationally a sample of M dwarfs with rotation axis having 
large inclination angle with respect to the line-of-sight (ideally per¬ 
pendicular) must be chosen and compared. Recent development 
of Zeeman-Doppler Imaging (ZDI) has provided us with magnetic 
maps of a number of stars. These maps can already be used t o tes t 
at least some dynam o scenarios. For insta nce, |Morin_gt_a 1. t 2()08a ! l . 
iMorin et~aD d2008bl) . iDonati et al.l d2008l) . and Morin et alT i 20 id) 
explored the surface magnetic field geometry and its evolution in 
a sample of M dwarfs consisting of partly and fully convective ob¬ 
jects. Only circularly polarized Stokes V spectra were used in this 
analysis. Therefore, the magnetic maps presented in these studies 
reflect mostly the large-scale magnetic structures that do not cancel 
out in polarized light (if a sufficient number of observed time-series 
is available to cover the full rotation cycle of the star uniformly). 

The main results of the ZDI studies are the existence of a 
switch in the magnetic field geometry between partly and fully 
convective M dwarfs: partly convective stars tend to host non¬ 
axisymmetric fields with dominant toroidal components, whereas 
axisymmetric dipole-like fields are found in fully convective stars. 
The distribution of field topologies, however, was found to be 
even more complicated. In particular, objects with weak toroidal 
fields were found among rapidly rotating fully convective stars. 
This pointed towards a possible dipolar breakdown so the stars 
with similar stellar parameters can host either domina nt polo idal 
or toroidal configurations (see Fig. 15 in iMorin et alj J20ld) ). A 
possible explanation for these observational findings was discussed 
in lGastine et al. I J 20131) who interpreted the magnetism of late M 
dwarfs as a result of a dynamo bistability and the choice of the final 
dynamo state in the models discussed above depends on initial con¬ 
ditions and does not change with time. We note here that, because 
our dynamo models are linear, they cannot provide evidence for 
bistability: the choice of the final dynamo state is fixed by the ini¬ 
tial conditions. It is important to stress that no signature of a global 
magnetic topology change in individual stars has been seen. It may 
be that much longer monitoring is required to observe such changes 
(which, in fact, would be a crucial test for the bistability scenario). 
Note however that the activity cycles in M d warfs are expect ed to 
be much longer than that of the Sun. lKitchatinov et al. I ( 120141) sug¬ 
gested a model for the observed change of state that depends on the 
underlying activity cycle. 

We note here some furth er points c oncerning the bi s tabilit y 
problem. iGastine et al. I d2013l) (see also ISchrinner et al. I J20 1 2| |) 
reported dynamo bistability at low Rossby Ro number from a 3D 
model, i.e. the dynamo action either produces a strong stable dipo¬ 
lar field without any significant time variation (ar 2 -dynamo), or a 
much weaker field which may possibly oscillate and be described 
by an a 2 D-dynamo. The choice of the final stable attractor depends 
on the initial conditions, or the prior physical conditions in the 
stars in their early stages after their formation. For slower rotators 
(Ro > 0.1), a dipolar breakdown is reported as strong dipole field 
become less likely when rotational influence on the convective flow 
becomes secondary. In principle, our illustrative model includes the 
possibility of a steady axisymmetric dynamo (dotted line in Fig. 1) 
and we recognize that further investigation may result in identifi¬ 
cation of this possibility with the actual field configurations in M 
dwarfs. 
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We can now try to compare our model predictions with the 
available ZDI maps. We first note that almost all stars with dom¬ 
inant non-axisymmetric toroidal components from the final sam¬ 
ple of lMorin et alj 1201(1) demonstrate similarities with the dipole¬ 
like, non-axisymmetric A1 mode (see Fig. 0. However, the ob¬ 
served orientation of the meridional component of the magnetic 
field does not agree with that predicted. For example, the radial 
field component of EV Lac is consistent with that of the A1 mode 
(see [Morin et alj|2008bl ). but the B^ component does not seem to 
show the expected surface distribution. If we ignore the B$ compo¬ 
nent and concentrate only on B r , then the signatures of mode A1 
(which changes in sign along longitude) are found among both the 
groups of partly and fully convective stars, such as, say, GJ 1245 B, 
GJ 1156, DT Vir, GJ 182 and some others. Note that the A1 mode 
is only a plausible explanation of the observed picture of the mag¬ 
netic field distribution in the framework of the dynamo modes that 
we consider in our study. Of course, this in no way proves that the 
actual magnetic field pattern in these stars can be described by such 
a simple model. Also note that our dynamo model may be not fully 
adequate especially at the stellar surface and that we show toroidal 
magnetic field distributions slightly below it. (Remember that our 
boundary condition is of zero toroidal field at the surface - see the 
short discussion in Sect. 2.) 

On the other hand, the mode S1 never seems to be observed. 
Many objects demonstrate a distribution of radial field, which 
is characterized by opposite polarities in Northern and Southern 
hemispheres at all longitudes, while a change in polarity would be 
expected. Thus this disagrees with predictions for both A1 and SI 
modes. Two stars, YZ CMi and DX Cnc, seem to show a possible 
transition between A1 and SI modes during subsequent observing 
epochs, but this results is not well constrained by the ZDI maps, 
and also the distribution of B^ does not agree with that of any of 
the predicted modes. 

The detection of AO and SO modes is more difficult because 
the existing ZDI maps better represent the large-scale magnetic 
structures and not the small-scale structures that might be asso¬ 
ciated with stellar spots. However, we note that objects with de¬ 
tected non-axisymmetric toroidal configurations often show local¬ 
ized magnetic structures that could be characterized as spots. Stars 
such as DS Leo, DT Vir, GJ 182, GJ 1156, and GJ 1245 show spotty 
patterns of the distribution of B, on their surfaces. These spotty 
structures evolve on time scales of years, and thus can be used to 
construct butterfly diagrams. Unfortunately, similarly to the cases 
discussed previously, the observed distribution of B^ does not agree 
with those predicted. In some stars we see the appearance of sin¬ 
gle large scale spots and not a group of spots with different polari¬ 
ties distributed along a given latitude, while in others (e.g. DT Vir, 
DS Leo, and OT Ser) very complex surface structures are detected 
and these stars are certainly good candidates for future monitoring 
and the detection of A0/S0 modes. In this context it looks attractive 
to suppose that the spots at M dwarfs are caused by radial magnetic 
field rather by eruptions of subsurface toroidal magnetic field, as is 
the case on the Sun. 

Another distinct feature of the observed magnetic field distri¬ 
butions is the presence of polar spots in the majority of fully con¬ 
vective objects. These spots can have positive (DS Leo, EQ Peg A) 
or negative (CE Boo, AD Leo, YZ CMi, WX UMa, GJ 51) signs 
of B r . It is unlikely that we observe a group of spots because only 
a single polarity is detected. Therefore, these magnetic polar spots 
look simply like poles of the large-scale dipolar field. Note that 
these polar spots also change their shapes, as detected for e.g., 
YZ CMi. This possibly implies that also other stars with polar 


spots should show the same evolution and their non-detection is 
connected with too short monitoring time. 

It is interesting to note that the fo rmati on of polar spots has 
recently been modelled bv lYadav et al.l J 201 5l) who found that, un¬ 
der certain conditions, their distributed dynamo models can spon¬ 
taneously generate large-scale dark spots at high stellar latitudes. 
Thus, these models represent an alternative explanation of the ob¬ 
served magnetic field geometries in rapidly rotating low mass stars. 

We may conclude at this point that the signatures of AO and SO 
modes, i.e. the presence of equatorial spots, can indeed be found in 
M dwarfs with strong toroidal fields, but it is hard to see stable pat¬ 
terns in their evolution. This is plausibly because the observations 
of these stars do not extend over long enough time intervals. 

It follows from the above considerations that it is not possi¬ 
ble to distinguish unambiguously between different dynamo modes 
from the available ZDI maps. The main limitation is the absence of 
long-time monitoring (i.e over several years). On the other hand, it 
is already possible to see a clear distinction between, say, A1 and 
SI modes in some stars, at least from the analysis of the radial field 
component. Additional monitoring is strongly needed. 

In order to put tighter constraints on the possible dynamo 
modes, observations of stars with large inclinati on angles (60° - 
70°) are needed. Among stars from lMorin et alj ( I2010I) . there are 
many that already meet this requirement, but none of them has an 
inclination larger than 70°. Nevertheless, the construction of the 
sample of stars with sufficiently large inclination angles does not 
seem to be problematic. The only difficulty is that inclination angles 
are often derived from ZDI itself so it is impossible to know these 
angles before actual phase-resolved observations. If rotational peri¬ 
ods are known (say, from analysis of observed light curves) then in¬ 
clinations can be derived from spectroscopically known v sin i val¬ 
ues. In many cases, however, spectroscopic usini values are not 
accurately constrained which thus results in large errors in incli¬ 
nation. Interferometry is an alternative technique that may provide 
model independent estimates of inclination angles, but its abilities 
(in most cases) are not yet sufficient to resolve the surfaces of small 
and faint M stars. 

Despite the detected changes in the magnetic field topology 
from Stokes V spectra, no strong evidence of the same changes 
have been n otice d from the analysis of unpolarized light, as dis¬ 
cussed in iShulvak et alj d2014h . The authors explored the distribu¬ 
tion of magnetic filling factors on surfaces of four M dwarf stars 
(from the very magnetically sensitive lines of the FeH molecule) 
and found similar distributions in all of them. Note that no infor¬ 
mation about the magnetic field orientation could be derived from 
this study. The only trend the authors could detect was in the dis¬ 
tribution of filling factors and the strength of the surface magnetic 
field, which seem to depend on the rotation rate. However this is 
still inconclusive because of the small sample of stars used. 

The only way to distinguish unambiguously between dynamo 
modes is by monitoring individual M dwarfs in polarized light 
by all possible means, in order to construct time-series of ZDI 
maps. By tracking the evolution of magnetic (and possible tem¬ 
perature) structures the corresponding butterfly diagrams could be 
constrained and comparison made directly with model predictions. 
Fig.[2]tells us that the monitoring should be quite long, over about 
30 years (to be compared with the 11-year solar cycle); 15-year 
monitoring would be sufficient for preliminary conclusions. 

Normally, ZDI techniques allow constraint of a few spherical 
harmonics of the surface magnetic field, with a typical maximum 
spherical harmonic degree of about 5-10. Thus, if the dipole has a 
typical oscillatory cycle of 60 years then it is very plausible that 
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higher degree modes can vary on shorter timescales. For instance, 
in geodynamo (a 2 ) models, the ti mescales of sec ular variations of 
the magnetic field obey a l/< law dchristensen et al.ti2012 ). High- 
quality observations and time-tracking over 5-10 years would then 
be sufficient to see a notable secular variation of the magnetic field 
in the higher spherical harmonic degrees. It looks thus promising in 
this respect to observe rapidly rotating stars such as, say, V374 Peg 
(P= 0.44d, iDonati et al. 1120061 : iMorin et at1l2008al) . Although no 
change in the magnetic field topology was reported between the 
two years of observations (2005 and 2006), continuing monitoring 
of this star is strongly recommended because the possible cycle is 
expected to be short. 


On the other hand, a significant amount of information con¬ 
cerning magnetic configurations can be obtained using DI only, cf. 
that the solar magnetic cycle was isolated by Schwabe in the XIXth 
century before Zeeman splitting in sunspots was observed. DI tech¬ 
niques have existed for about 30 years but the time resolution of 
the published maps of M dwarfs at best only span several years. 
Note that the length of activity cycles is expected to scale with ro¬ 
tation rate. This_issue was seems to have been first discussed by 
iNoves et al. I dl984l) . The tendency is to some extent confirmed b y 
observations in the Wilson program (Saar & Br andenburg II 199^1 . 
Observing the most rapidly rotating stars may be a relatively easy 
task. For instance, the rapidly rotating YZ CMi, DX Cnc, GJ 1156, 
GJ 51, etc. are potentially good targets with which to begin. In 
addition, the strength of the surface magnetic field clearly corre¬ 
lates with the rotation rate: it grows steadily with decrease of rota¬ 
tion period until a saturation limit is reached (Reiners et all l2009h . 
Importan tly, as was shown inlPizzolato et al] d2003l) and more re¬ 
cently in iReiners et al. !B, the magnetic field strength in the 
non-saturated regime is a function only of rotation rate and does 
not depend on any other stellar parameters (such as mass or ra¬ 
dius), while the saturation limit depends on the bolometric lumi¬ 
nosity and thus differs from star to star. Therefore, studying rapid 
rotators could produce interesting constraints on dynamo action in 
M stars. 


So far all ZDI maps have been obtained by using ob¬ 
servations in Stokes V only. In order to resolve more sur¬ 
face detail, mapping of both linear and circular polarization is 
needed. This imposes very strong observational constraints be¬ 
cause, e.g., long integration times would be needed to collect 
sufficient signal in linearly polarized light. High spectral reso¬ 
lution is also an essential requirement. The instruments avail¬ 
able for M dwarf research are ESPaDOnS@CFHT (3.6 m, K = 
65 000, 370-1005 nm), NARVAL @TBL (2 m, R = 65 000, 370- 


1005 nm) dDonati et alj Il997l). HARPSpol@ESO (3.6 m R ~ 
110000 , 378-691 nm) dSnik et al] I20T1I : IPiskunov et alj 1201 lh . 
and HiVIS@AEOS (3.7 m, R « 50000, 500-1000 nm; R « 
33 000, 1000-2500 nm) dThomton et alj l2003l : lHarrington et al.1 
120061: lHarrington & Kuhnl l2008l) ~ as well as future instruments 
such as PEPSI@ LBT (2 x 8.4 m. R = 120000, 383-907 nm) 
dllvin et alj 1201 lh. SPIRou@C FHT (3.6 m, R » 70000, 980- 
2350 nm) lArtieau et al.l 1201 ll), and CRIRES @VLT (8 m, R = 
100000, 1000-5300 nm) dKaeufl et aid 120041) . Observing at in¬ 
frared wavelengths with high spectral resolution is superior to vi¬ 
sual observations because of the cool temperatures of M dwarfs, 
and that the scaling of Zeeman splitting is proportional to A 2 . Thus 
instruments operating in the near-infrared would contribute greatly 
to studying the evolution of magnetic topologies in M dwarfs. 


5 DISCUSSION AND CONCLUSIONS 

Based on a relatively simple and standard mean field dynamo 
model for M dwarfs we have presented four dynamo generated 
magnetic configurations which can be excited in the framework 
of the model. In this model, the parameter which determines the 
excited dynamo mode is the magnetic Prandtl number P m . Note 
that we identify the excited mode with the mode with the lowest 
marginal dynamo number. We appreciate that for more supercrit¬ 
ical non-linear models this simple idea may become inadequate. 
However, we intend to begin our analysis with simple models and 
improve them if observations will so require. 

We started from the conventional expectation that observa¬ 
tions and dynamo theory are still quite remote from each other, 
and comparison is a very uncertain and problematic undertaking. 
After making the comparison it appears plausible that the situa¬ 
tion may not be quite as difficult as expected and some preliminary 
identifications may be obtained immediately. We see that at least 
some of the configurations look closer to observations than oth¬ 
ers. In particular, signatures of the SI configuration appear never 
to have been observed, while prospects for identification of the Al 
configuration perhaps look more promising. The surface distribu¬ 
tion of radial magnetic field looks rather similar to the predictions 
of dynamo theory while the distribution of toroidal components in 
dynamo models appears quite different to the observational data. 

Note that systematic investigation of magnetic activity of M 
dwarfs can be useful for the understanding of stellar hydrodynam¬ 
ics. Indeed, Eq. 0 suggests that dynamo action becomes stronger 
for stars with larger stellar radius, if other governing parameters are 
unchanged. In practice however turbulent diffusivity increases with 
stellar radius so the stellar hydrodynamical model predicts that dy¬ 
namo action becomes weaker for more massive stars. On the other 
hand, Fig. |T| tells us that preferred magnetic fields configurations 
are specific for weaker and stronger dynamo action, as given by 
the C a parameter. Confronting these expectations with future ob¬ 
servations we may hope to deduce the actual scaling of turbulent 
diffusivity with stellar radius. 

In principle, we could start fitting dynamo models in order to 
reproduce an improved phenomenology for However our feel¬ 
ing is that, taken overall, the results on the topic are still not stable 
enough, and it would be better to discuss the available and forth¬ 
coming results in the framework of the strategy suggested here, as 
well as bearing in mind the possibilities discussed above when an¬ 
alyzing any forthcoming and improved dynamo modelling. In any 
case, it looks plausible that the joint effort of observers and dynamo 
modellers will be able to clarify the problem in the foreseeable fu¬ 
ture. 

We recall the brief overview of dynamo modelling given in 
Sect. |T| and appreciate that the dynamo model we have used can 
of course be criticized on a number of grounds. However it seems 
remarkable that all dynamo driven magnetic field configurations 
for M dwarfs obtained in its framework give a magnetic field that 
is concentrated at quite high latitudes, i.e. they predict magnetic 
stellar activity in regions that are far from the stellar equator. Based 
on solar experiences, this result is quite unexpected and deviates 
substantially from expectations. On the other hand, the formation of 
high-latitude spots in rapidly rotating fully-convective stars is also 
predicted bv self-consistent global dynamo models presented by 
lYadav et al.l d2015l) . This may indicate that the dynamo mechanism 
in cool rapidly rotating stars is quite different from the solar case. 

We finally stress again that our intention has been to explore a 
possible methodology for investigation of the phenomenology of M 
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dwarf magnetism. We have illustrated this by reference to a particu¬ 
lar dynamo model, but as more sophisticated and definitive models 
for the magnetic fields are developed, the principles will remain 
relevant. 
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